Draft version February 2, 2008 

Preprint typeset using LATpjX style cmulateapj v. 11/12/01 



TEV EMISSION FROM THE GALACTIC CENTER BLACK-HOLE PLERION 

Armen Atoyan 1 and Charles D. Dermer 2 

1 CRM, Univcrsitc dc Montreal, CP. 6128, Montreal, Canada H3C 3J7; atoyan@crm.umontreal.ca 
2 EO.Hulburt Center for Space Research, Code 7653, Naval Research Laboratory, 
Washington, DC 20375-5352; dermer@gamma.nrl.navy.mil 
Draft version February 2, 2008 

ABSTRACT 

The HESS collaboration recently reported highly significant detection of TeV 7-rays coincident with 
Sgr A* (Aharonian et al. 2004). In the context of other Galactic Center (GC) observations, this points 
to the following scenario: In the extreme advection-dominated accretion flow (ADAF) regime of the 
GC black hole (BH), synchrotron radio/sub-mm emission of ~ 100 MeV electrons emanates from an 
inefficiently radiating turbulent magnetized corona within 20 Rs (Schwarzschild radii) of the GCBH. 
These electrons are accelerated through second-order Fermi processes by MHD turbulence, as suggested 
by Liu et al. (2004). Closer to the innermost stable orbit of the ADAF, instabilities and shocks within 
the flow inject power-law electrons through first-order Fermi acceleration to make synchrotron X-ray 
flares observed with Chandra, XMM, and INTEGRAL. A subrelativistic MHD wind subtending a ~ I sr 
cone with power > 10 37 ergs -1 is driven by the ADAF from the vicinity of the GCBH. As in pulsar 
powered plcrions, electrons are accelerated at the wind termination shock, at > I0 16 5 cm from the 
GCBH, and Compton-scatter the ADAF and the far infra-red (FIR) dust radiation to TeV energies. 
The synchrotron radiation of these electrons forms the quiescent X-ray source resolved by Chandra. The 
radio counterpart of this TeV/X-ray plerion, formed when the injected electrons cool on timescales of 
> 10 4 yrs, could explain the origin of nonthermal radio emission in the pc-scale bar of the radio nebula 
Sgr A West. 

Subject headings: acceleration of particles — black hole physics — Galaxy: center — plasmas — radiation 
processes: nonthermal 



1. INTRODUCTION AND OBSERVATIONS 

HESS measurements of TeV fluxes from Sgr A* (Aha- 
ronian et al. 2004) confirm, with high significance, earlier 
reports by the Whipple (Kosack et al. 2004) and CAN- 
GAROO (Tsuchiya et al. 2004) collaborations (albeit at 
higher flux levels than HESS) about TeV emission from 
the Galactic Center. This result adds important new in- 
formation about the massive (Mgc ~ 4 x 10 M©; Chez 
et al. (2003)) BH at the center of the Galaxy. When com- 
bined with the wealth of observational data in the radio 
to X-ray wavebands (for review, see Melia & Falcke 2001), 
the TeV 7-ray data point to a scenario involving an inner, 
inefficiently radiating magnetized corona, i.e. the ADAF, 
and a subrelativistic outflow of particles and field from the 
ADAF. As in pulsar plerions like Crab Nebula (see Ken- 
nel & Coroniti (1984)), the wind terminates at a shock, at 
R > 3 x 10 16 cm from the GCBH, to create a black-hole 
powered 'plerion' where the TeV radiation is produced. 

The VLBI size of the radio source at the location of 
Sgr A* is 110(±60) ^as at v = 215 GHz and 190(±30) //as 
at v = 87 GHz (Krichbaum et al. 1998). Using VLBA, 
Bower et al. (2004) resolve Sgr A* at 42 GHz to have an 
angular size of 240(±10) /z, or R rad ~ 24 R s for M GC 
- x 1O 6 M black hole and a GC distance d = 8kpc, and 
find that the size increases with wavelength. The radio 
emission of Sgr A* is rather stable, with total luminosity 
£36 = L rac i/10 36 ergs -1 ~ 1 and spectral flux S v oc u 3 
that sharply falls off above v$ ~ I0 12 Hz (Zylka ct al. 1995; 
Falcke et al. 1998). 



The near-infrared (NIR) flares show factors-of-2 flux 
variations on timescales of ~ 20 min and reach powers of 
w 10 35 ergs -1 (Genzel et al. 2003). The spectra rise in vS v 
diagrams, and show powers w 10 35 ergs _1 , similar with 
the (nonsimultaneously) detected X-ray flares (Baganoff 
et al. 2001; Porquet et al. 2003). The 200s X-ray variabil- 
ity timescale limits the engine size to r — R/Rs < 5. The 
unabsorbed quiescent flux detected by Chandra observa- 
tory from Sgr A* is L x . gui (2-10keV) « 2.4 x 10 33 ergs~\ 
and has angular extent #fwhm ~ 1-4" (Baganoff et al. 
2003), corresponding to source size Rx ~ 8.3 x 10 16 cm. 

The TeV spectral fluxes in two HESS campaigns in 
June/ July and July/ August of 2003 are almost exactly the 
same, extending with photon index 2. 1-2.3 from 160 GeV 
to 10 TeV, with total power L 7 « 10 35 ergs _1 . There is 
no sign of flux variation within the statistics available in 
each of these sets (Aharonian et al. 2004). Given the high 
statistical significance of the signal in both data sets, these 
results suggest a quiescent or quasi-stationary TeV source 
coincident within < 1.5' of Sgr A*. 

We propose a model for Sgr A* where the quasi- 
stationary radio and flaring X-ray/NIR emissions are syn- 
chrotron radiation from the ADAF. The MHD wind from 
the ADAF powers a black-hole plerion where the quiescent 
X-ray and TeV emissions are produced by electrons accel- 
erated at the wind shock. The model explains the existing 
multi-wavelength data, and makes predictions that can be 
checked by various X-ray and gamma-ray detectors. 

2. PARAMETERS OF THE ADAF 
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The dimcnsionless mass-accretion rate (or Eddington ra- 
tio) m = tibhMc 2 / LEdd is denned in units of the Edding- 
ton luminosity LEdd, where t/bh — 0.1 is the maximum 
efficiency of gravitational-to-photon energy conversion by 
the BH. In the ADAF model (Narayan & Yi 1995; Esin et 
al. 1998), thermal radiation of the accretion plasma at the 
level Lt = m 2 LEdd/rn* is predicted when rh < rh* ~ 0.1. 
The rate rh ~ 1.5 x 10~ 5 is obtained if one formally equates 
Lt with the observed L ra d — 10 36 ergs _1 . 

The compact radio emission cannot originate from an 
optically-thick accretion disk which would only be allowed 
at R ^> Rrad, or from a much hotter ADAF plasma at 
smaller scales, but is compatible with a synchrotron origin 
(Liu, Petrosian, & Melia 2004). The thermal output of 
the ADAF would be peaked mostly in the hard X-ray/soft 
gamma-ray domain, where the quiescent luminosity is at 
least 2 orders of magnitude smaller than L ra d- This results 
in the accretion rate estimate m_6 = m/10~ 6 < 10~ 6 con- 
sistent with ADAF model and Sgr A* observations. Even 
for rh this small, < 10 39 ergs _1 of accretion power is still 
cither advected into the BH or escapes as an outflow. 

Adopting GCBH mass 3 x 10 6 M , the ADAF magnetic 
field at radii r = R/R s is (Narayan & Yi 1995) 

B(r) « ^Qa-^X^m^r-^ G , (1) 

where A = 1 — (3 is the ratio of the magnetic to the total 
pressure in the accreting gas and a is the viscosity param- 
eter. Using a ~ 0.1 and A ~ 0.15 (Quataert & Narayan 
1999), the magnetic field B ~ 10 G in the ADAF at r < 20 
is found. For these fields, the Lorentz factor of electrons 
producing synchrotron radiation peaked at i/q > 10 12 Hz 
is 70 ~ 200, and the ADAF is still transparent to self- 
absorption for >GHz radiation. 

A stochastic acceleration model for the radio-mm emis- 
sion from Sgr A* (Liu, Petrosian, & Melia 2004) can 
be quantified using the second-order Fermi acceleration 
theory of Dermer, Miller, & Li (1996). The mag- 
netic field is defined through the relation B 2 /8n = 

sbVbhMc 2 /(4tt.R 2 c), implying B « 30eg /2 L^ 4 G for a 
region of size 20i?s- Equating the acceleration rate of 
electrons by whistler turbulence with the synchrotron loss 
rate gives a characteristic electron Lorentz factor 70 = 

5 (C £ s) 1 ^ 3 ^36 6 ' r T 11 ^ 18 , assuming a Kolmogorov spectrum 
for the turbulence. Here £ is the fractional whistler tur- 
bulent energy density compared to B 2 /8tt, and tt is the 
Thomson depth of the ADAF, assumed to be composed of 
e-p plasma. Values of B = 10 G and 70 = 200 are implied 
for the ADAF when C « 0.1, e B w 0.1 and t t « 2 x 10~ 4 . 

3. BLACK HOLE PLERION 

A Compton (i.e., leptonic) origin for non- variable TeV 
flux detected with HESS implies that the synchrotron 
counterpart of that flux is the quiescent component of the 
X-ray flux. The magnetic-field energy density in the source 
should then be ti^ = B 2 /8n < 0.1 u ra d, where 

urad = = 2.65 x 10- 8 ||crgcm- 3 (2) 

is the energy density of the target photons at R = 
i?i6l0 16 cm. Inside the radio source, at r < 20, the 
radiation energy density saturates at the level u ra d — 
L ra d/2TrcR 2 ad . Thus, the origin of TeV flux at those small 



distances would require magnetic fields B ~ 0.1 G, far be- 
low the ADAF model values. 

In the ADIOS ( "advection-dominated inflow-outflow so- 
lution" Blandford & Begelman (1999)) extension of the 
ADAF model, the outflow may carry a significant frac- 
tion of the generated kinetic energy of the plasma from 
the BH vicinity to large distances. We assume a wind of 
magnetized accretion plasmas with total power L win d = 
L w .37 10 37 ergs _1 propagating in a two-sided cone with 
opening angle Q at speed v w < c; for calculations we use 
v w — c/2. This wind terminates at a subrelativistic shock 
at Rshock < Rx from the GCBH. The shock accelerates 
electrons to energies E max = 7 maa; m e c 2 3> 10 TeV to pro- 
duce the observed TeV 7-rays by Compton upscattcring 
of sub-mm (v a ~ 10 12 Hz) photons from the ADAF. An- 
other important photon target is the emission from the 
cold (~ 100 K) dust ring of Sgr A West, with total lumi- 
nosity ~ 5 x 10 6 Lq in the central arcmin (few pc) region 
around Sgr A* (Bccklin et al. 1982). This radiation has 
a spectral peak at ~ 10 13 Hz and estimated energy den- 
sity ~ 2.4 x 10 3 eVcm -3 . It becomes the main contribu- 
tor to the < 100 GeV Compton flux from the plerion at 
R > 10 17 cm from Sgr A*. Synchrotron radiation from 
these same multi-TeV electrons at R < 10 17 cm produces 
the quiescent X-ray flux from Sgr A*. 

Chandra observations (Baganoff et al. 2003) reveal dif- 
fuse thermal X-rays with luminosity 2.4 x lO^ergs" 1 
from the central parsec region of Sgr A*, produced by 
kT ~ 1.3 keV plasma with n e = 26/rjy 2 cm^ 3, density, 
where r\f = ?7_i/10 _1 is the volume filling factor. Equat- 
ing the gas pressure n e kT in that region with the energy 
density of the wind gives the shock distance 

R shock ~ 3.1 x 10 16 L^Sl- 1 ' 2 ^ cm. (3) 

The TeV luminosity L 1 ~ 10 35 ergs" 1 implies a total elec- 
tron acceleration power L e > 10 36 ergs _1 . Assuming that 
the efficiency of electron acceleration L e / L W i n d < 30%, 
and letting the rest go to thermal plasma and nonthermal 
protons (unlike in pulsar winds where relativistic electrons 
dominate), we arrive at £^.37 > 0.3 . Note that Eq.(3) pre- 
dicts that the maximum power of the outflow sustainable 
at distances R s hock < Rx is L w ,^ < 10. 

The mean fluid speed Vf in the TeV plerion formed 
downstream of the shock can be derived from (Rx — 
Rshock) ^ Vf x t c (lx), where 

Mix) , 4x1 „ 7 ||(_m_)-' 2 s (4) 

is the Compton cooling time of electrons producing syn- 
chrotron X-rays with energy ex for the estimated mag- 
netic field Bpi er ~ 100 /xG in the GCBH plerion at scales 
R < Rx- Taking i? 16 = 5 for the median distance of the 
TeV plerion, we find vj > 500km/s at those distances. 

Since relativistic electrons of all energies convect with 
the fluid at the same speed, the TeV plerion model pre- 
dicts shrinking of the quiescent X-ray source size at higher 
energies. Note also that the opening angle of the outflow 
should be ft ~ 1 rather than ft <C 1, since otherwise R s hock 
would exceed Rx- This indicates that the energy outflow 
is in the form of a wind rather than well-collimated jet. 
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4. SPECTRAL MODELING 

Fig. 1 shows numerical calculations of the quiescent ra- 
diation components of Sgr A* . These include synchrotron 
and SSC emission from the magnetized ADAF within 
r ra d ~ 20 of the GCBH {dashed curves), X-ray and TeV 
emission from the compact plerion at scales R ~ (3 - 
10) x 10 16 cm (solid curves), and the emission from the 
larger plerion (dot- dashed curves) inflated to pc scales in 
the process of convective (and possibly also diffusive) prop- 
agation of the accelerated electrons on timescales < 10 4 yr. 
Magnetic field in the latter, B^ = 140 [iG, is assumed 
higher than B\ = 90 /xG in the TeV plcrion. This is in 
agreement with the expected increase of the magnetic field 
downstream of the shock in typical pulsar plerions (Ken- 
nel & Coroniti 1984), and is explained by deceleration and 
compression of the plasma. But it is also possible that 
£?2 < B\, in which case a smaller synchrotron flux from 
the pc-scale plerion than shown in Fig. 1 is predicted . 

The inset in Fig. 1 shows the nonthermal electron distri- 
butions. Besides synchrotron and Compton losses, calcu- 
lations also include Coulomb and brcmsstrahlung losses in 
a medium with n gas = 10 3 cm~ 3 characteristic for central 
pc regions of Sgr A*. Electrons with a power-law injec- 
tion index a p i er — 2.2 and an exponential cutoff above 
E m ax = 50 TeV are injected into 2-sidcd nebula with total 
power L e = 6.5 x 10 36 ergs _1 . These electrons radiatively 
cool downstream of the shock and propagate to distances 
> 10 17 cm on the timescale t esc — 50 yr, during which most 
of the injected energy of the mufti- TeV electrons (produc- 
ing TeV and X-ray fluxes observed) has been already lost. 

The radio emission is produced at r < 20 in B « 10 G 
field by a relativistic Maxwellian electron distribution 
oc -y 2 e~ 7//70 with 70 = 200 resulting from a balance be- 
tween second-order energy gains and synchrotron losses 
(Schlickciscr 1985). The stochastic acceleration power of 
radio electrons is therefore about equal to the observed 
radio luminosity L rac i- 

The SSC emission component of these electrons, shown 
in Fig. 1 by the dashed curve, falls in the sub-keV region 
below the level of the quiescent X-ray flux. However, if 
magnetic fields < 5 G are assumed, this SSC contribution 
would rise and also move to 1 keV, and could contribute 
to radiation at that energies. It cannot, however, explain 
the entire quiescent flux extending to 10 keV. 

The powerful X-ray and NIR flares observed from 
Sgr A* are explained in our model by the onset of in- 
stabilities of the accretion flow at distances of only few 
Rs- They result in strong shocks in the accretion flow 
and effective acceleration of particles advected with the 
flow. The accelerated particles are then easily taken out 
of this region in the wind/jet outflow. Synchrotron emis- 
sion of electrons accelerated to 7 > 10 6 explains the X- 
ray flares with very short variability scales. Furthermore, 
while propagating through the first few tens of Rs, these 
electrons have sufficient time, <~ 10 3 s, to cool in the high 
B fields there down to 7 < 3 x 10 3 to produce flares in the 
NIR domain. Self-absorbed flares at < 100 GHz detected 
on < 1 day timescales after the X-ray flares (Zhao et al. 
2004) could be explained by the radiation from these same 
electrons at later stages/larger distances of the outburst in 
the 'expanding source' scenario. 

In Fig. 2 we show the flare fluxes expected in this model. 



After shock-acceleration close to the BH and injection into 
collimated wind outflow with speed c/2, electrons prop- 
agate through the inner r ~10-30 region on timescale 
t esc — 1200 s, after which the radiation losses drop because 
of wind expansion and decline of B. The solid, dashed, 
and dot-dashed curves show synchrotron and Compton 
fluxes produced at times 200 s, 1200 s and 1 hr after the 
onset of the flare. The electron injection time profile is 
L e .fiare(t) = L /(l + t/t ) 2 , with t = 720 s. The mean 
magnetic field in ADAF outflow could be enhanced at 
the flaring state, so we take B = 25 G. The Compton 
fluxes shown (thin curves) clearly demonstrate that no de- 
tectable TeV flares should be expected during powerful 
X-ray flares. This is in agreement with the non-detection 
of TeV flux variations during many days of observation 
with the HESS telescopes, whereas the X-ray flares occur 
with frequency of ~ 1 per day (Baganoff et al. 2003). 

5. SUMMARY AND CONCLUSIONS 

A model consisting of magnetized coronal ADAF (syn- 
chrotron) radio emission within ss 20 Rs of the GCBH, 
and of a black-hole plerion powered by a wind from the 
ADAF resolves many puzzling observations of Sgr A* . X- 
ray flares are synchrotron radiation of electrons acceler- 
ated through first-order Fermi process by shocks within a 
few Rs of the GCBH. Electrons accelerated to 7 > 10 8 at 
the wind termination shock at > 3 x 10 16 cm creates the 
GCBH plcrion. Its synchrotron X-ray emission has been 
resolved as an ~ 1.4" source with Chandra (Baganoff et al. 
2003). The mufti- TeV electrons Compton scatter the radio 
photons from Sgr A* and FIR photons from the dust ring 
of Sgr A West to produce the TeV emission. TeV emission 
from the GCBH plcrion is nearly stationary because the 
cooling time of TeV electrons is ~ 100 yrs. 

A "jet-ADAF" model for Sgr A* has been proposed by 
Falcke & Markoff (2000) and Yuan, Markoff, & Falcke 
(2002). Our model, though also based on energy out- 
flow from the ADAF, differs greatly. In particular, the 
quiescent and flaring X-ray components are produced at 
different sites. The origin of the X-ray and NIR flares 
are explained as synchrotron emission of electrons acceler- 
ated during episodes of instabilities very close to the BH, 
not as Compton radiation in the SSC scenario of the "jet- 
ADAF" model. Importantly, the TeV flux cannot be eas- 
ily explained if the observed X-rays were due to the SSC 
mechanism in the vicinity of the GCBH, as suggested by 
Falcke & Markoff (2000). 

Propagation of GeV electrons from the plerion on 
timescales of > f0 4 yr with speeds ~ 100kms _1 could sig- 
nificantly contribute to the radio synchrotron flux of Sgr 
A West. In particular, it could explain the suggested 
nonthermal origin of radiation in the bar of Sgr A West 
(Wright et al. 1987), formed in the outflow direction. 

We also predict that quasi-stationary Compton and 
brcmsstrahlung fluxes from the pc-scale plerion, coinci- 
dent with the central parts of Sgr A West, will be sig- 
nificantly detected and possibly resolved with GLAST at 
GeV energies. But the expected highly variable Comp- 
ton counterpart of the synchrotron X-ray flares from the 
GCBH vicinity is too weak to be detectable with GLAST 
or HESS. The synchrotron extension of the flare emission 
in the hard X-ray/soft 7-ray domain explains the ~ 40 min 



4 



Atoyan & Dermer 



episode of a profound increase of flux detected by INTE- 
GRAL from the direction of Sgr A* (Belanger et al. 2004). 

The absence of apparent TeV flux variations may sug- 
gest a proton origin for the TeV radiation (Aharonian 
& Neronov 2004). Indeed, proton and ion acceleration 
through first- and second-processes in the ADAF and 
through first-order acceleration at the wind termination 
shock could make cosmic rays to produce TeV emission 
through nuclear pp-interactions with n ~ 10 3 cm -3 dense 
gas on pc scales, and could form extended TeV emission 
possibly already detected with HESS (Aharonian et al. 
2004). The hadronic origin of the TeV radiation in the 
ADAF itself in the BH vicinity requires, however, an ex- 
tremely dense gas target or extremely large proton powers 
> HP'ergs- 1 . 

Finally, we note that the unidentified EGRET source 
3EG J1746-2851 towards the GC (Mayer-Hasselwander 



et al. 1998) is significantly displaced (Dingus & Hooper 
2002) from the direction of the GCBH, and is unlikely 
to be related to Sgr A*. It is probably emission from a 
young pulsar, though not with the "mouse" PSR J1747- 
2958 (McLaughlin & Cordes 2003). A young pulsar with 
7-ray properties like Vela but with apparent 7-ray power 
« 10 x larger could have been missed in pulsar surveys 
due to the large dispersion measure towards the GC. We 
suggest a deeper, higher radio-frequency and X-ray search 
at the refined location of 3EG J1746-2851. 

We thank I. Grenier, S. Markoff, R. Narayan, P. Ray, 
and R. Romani for discussions. Research of CD and visits 
of AA to the NRL High Energy Space Environment Branch 
are supported by GLAST Science Investigation No. DPR- 
S-1563-Y. The work of CD is supported by the Office of 
Naval Research. 
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Fig. 1. — Quasi-stationary/quiescent synchrotron (heavy curves) and Compton (thin curves) fluxes expected from Sgr A*. Dashed curves 
show the fluxes from the magnetized ADAF corona with B = 10 G and R ~ 20 Rs due to a relativistic Maxwellian distribution of electrons 
with mean Lorentz-factor 70 = 200. The solid curves show the fluxes from the BH plerion with B\ = 90 fiG at a distance R p i er ~ 5 X 10 16 cm; 
the calculations assume escape of the flow from the TeV plerion to R > 10 17 cm on timescale 50 yr. The dot-dashed curves show the fluxes 
formed during assumed t esc = 2 X 10 4 yr propagation time in the plerion at scales R ~ 0.03 — 3pc. The dotted curve shows the bremsstrahlung 
from the central few pes of Sgr A West with gas density n gaa = 10 3 cm -3 . The inset shows the electron energy distributions formed in those 
3 spatial scales. The hatched regions in the X-ray domain show the fluxes detected by Chandra (ax ~ 2.2 — 3 are plotted) and XMM 
(ax = 2.5 ± 0.3), and the TeV fluxes detected by HESS (see text for references). The radio data are from Zylka et al. (1995) and Falcke et 
al. (1998). 
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Fig. 2. — Synchrotron and Compton fluxes (thick and thin curves, respectively) expected during a flare at 200s (solid), 1200s (dashed), and 
3600 s (dot-dashed curves) after the onset of the injection of relativistic electrons with ct p i er = 2.2 and j m ax = 10 7 , initial injection power 
L e (0) = 6.3 X 10 36 ergs -1 , and characteristic decline time to = 720s (see text). Escape of electrons with the conical wind during t esc = 1200s 
from the central R < 20 Rs = 1.8 X 10 13 cm region of the ADAF is supposed. The dotted curve shows the expected bremsstrahlung 7 rays at 
t = 200 s, calculated for plasma density of 10 8 cm -3 in the ADAF. 



